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Abstract. Thin metal films are subject to the pressure caused by the zero point 
oscillations of the electromagnetic field, which depends upon the film optical properties 
and, in case of deposition onto a substrate, upon the substrate reflectivity. It has been 
suggested that this force may be relevant in determining the stability of deposited 
pseudomorphic films with respect to buckling or island formation. We present a 
detailed analysis of its behaviour as a function of the optical parameters and of the 
film thickness and we illustrate the conditions under which it may play some role. For 
free standing films it turns out that the film stabilization is basically due to the surface 
' stress, which largely overwhelms the vacuum force. For epilayers on metal substrate the 

CN| . vacuum force may be important, and we give stability diagrams and critical thicknesses 

for several cases, illustrating how the flat surface growth may depend upon the film 
parameters. The importance of including retardation effects into the theory for a 
realistic determination of the stability conditions is also discussed. 



PACS numbers: 42.50.Lc,03.70.+k,73.20.Mf 



1. Introduction 



Epitaxial deposition of thin metal films on various substrates is of crucial importance in 
several applications ranging from the design of systems for heterogeneous catalysis and 
electron emission to the technology of semiconductor devices and magnetic recording 
[H El El H]. In several cases a significant mismatch between the lattice parameters 
of the deposit and the substrate is present. If the mismatch is not too large, the 
overlayer grows in a way that the atoms of the deposit are in registry with the substrate 
atomic structure and the growth is pseudomorphic. As the film thickness increases, 
mechanical strain energy builds up in the deposited film, which eventually can cause 
instability unfavorable to uniform flat film growth: the formation of dislocations or the 
continuation of growth in the form of islands become energetically favorable. In several 
cases epitaxially deposited film grown layer by layer can undergo a transition to a growth 
mode of three dimensional dislocation-free islands, which form on top of a layer of a 
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certain critical thickness [5J EE El El E] • It turns out that the critical thickness is material 
dependent and decreases with increasing the lattice mismatch strain [T0| [TT]. 
The structure and morphology of the deposited films is often determined by kinetic 
constraints rather than by the minimization of the free energy. However it is possible to 
work in experimental conditions that allow the system to relax toward the equilibrium 
configuration. In such cases equilibrium considerations may be useful to understand 
the physical interactions that are responsible of the film evolution and to establish the 
conditions for the existence of a critical thickness and its value. 

One can discuss the stability issue using a continuum model [121 EE EEU EE EE EH]- The 
value of this approach lies in the fact that one can draw qualitative and some quantitative 
conclusions on general aspects of the growth without the need of a microscopical 
description of the atomistic mechanisms present in the process. Within this model 
the stability of epitaxial films results primarily from the competition between the 
surface energy and the elastic energy, caused by lattice mismatch, since when a flat 
film surface is modified into a wavy shape the surface energy increases, but the elastic 
energy decreases. Roughening of the surface and island formation can occur when it 
is energetically favorable to relax the elastic energy by increasing the surface area. It 
has been observed that since the surface energy tends to stabilize perturbations at low 
wavelength, while the elastic energy amplifies the perturbation at all the wavelength, 
any strained film, regardless of its thickness, should be unstable. This prediction does 
not match with the experiments which show the existence of a critical thickness beyond 
which the instability occurs. This suggests that some additional force has to be present 
that competes with elastic stress and stabilizes the film within a range of thickness of 
the order of few nanometers. Long range forces, in particular the vacuum fluctuation 
forces between the interface and the film free surface, have been considered as a possible 
source of stability in the system [101 [18]. Due to their peculiar size dependence, these 
forces are expected to play a role only at nanometric distances and when they cause a 
repulsion between the films boundaries. 

In this paper we discuss the conditions under which the vacuum fluctuation forces can 
contribute to deposited metal film stability. Our purpose is to identify, both for the film 
and for the substrate, the range of optical parameters, that can determine the existence 
of a finite critical thickness below which the flat surface growth takes place. To this 
end we need a realistic description of the metal optical properties. Recently we have 
reported on the results of a study of vacuum electromagnetic forces on thin metal films 
based on the plasma model [19] . We have shown that the force may change in sign 
and intensity depending upon the film size and the nature of the substrate. However 
the plasma model does not provide an accurate description of the real metal dielectric 
properties, since it neglects both the intraband absorption of the free electrons and the 
effects on the interband transitions in the relevant frequency range. 
A more accurate description is give by the Drude model 
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which includes the parameters tt p , the plasma frequency, and u T the relaxation 
frequency. For uj t = 0, corresponding to infinite relaxation time, we recover the plasma 
model. While in the plasma model Q p is fixed by the valence electron density, the 
Drude model allows to extract the parameters from the experimental data, thus taking 
into account both the real electron density and the effective mass of the electron in the 
metal under consideration. Recently a careful study of the fitting procedure for Au films, 
focused on the frequency range that is relevant in the determination of the forces due 
to vacuum electromagnetic field fluctuations, has shown that the Drude model can be 
useful to reproduce optical experimental data for films prepared in various experimental 
conditions [20], thus allowing for a more accurate evaluation of the force. 
We start our analysis with a study of the force on isolated metal films. A free 
standing metal film is subject to a negative pressure i.e. an attractive force 
per unit area F between the boundaries caused by zero point oscillations of the 
quantized electromagnetic field. We have performed calculations of the vacuum induced 
electromagnetic force at T = 0°K for free standing films of different thickness using the 
dielectric function of equation ([!]) and analyzed its dependence upon the optical film 
parameters. Here we report on the results and on related theoretical aspects concerning 
the role of this force in the film stability. The study of the forces on deposited films 
has been carried out in two steps. First we have considered the force on a metal 
film deposited onto a perfectly reflecting substrate. In this case the force between 
the film boundaries, i.e. the substrate-film interface and the film surface, is always 
repulsive, independently of the film parameters, and the discussion on the stability 
issue, carried out with reference to standard values of the stress concentration and the 
surface energy, is particularly simple. We identify the range of parameters that permit 
the existence of a critical thickness in this ideal case. Turning to real metal substrates 
we find, in agreement with previous work based on a plasma model description [19], 
that the force can be repulsive or attractive, depending upon the the difference in the 
optical parameters of the substrate and the film. An attractive force contributes to 
the destabilization of the flat surface growth and favors the island formation. In case 
of a repulsive force we provide examples of stability diagrams, which give the range 
of parameters that allow the system to remain stable, and we plot the corresponding 
critical thicknesses. Our purpose is not to derive results for a specific system, rather we 
illustrate the conditions under which this force may be important. 

2. Free standing metallic film 

The expression of the vacuum fluctuations energy at T = 0°K in a configuration with 
two semi-infinite media of dielectric functions ei(u>) and 62(0;) separated by a film of 
thickness d and dielectric function 63(0;) (see fig{T]) is given by the Lifshitz formula [21J. 

hi? f°° f°° r 1 
E ^ = ^J 1 pdp l t\HQTM{it)] + ln{Q TE {it)))dt (2) 
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Figure 1. Notation for three layers system 



where: 



(eiK 3 + e 3 Ki)(e 2 K 3 + e 3 K 2 ) 

^-*- £;g£;g «- w w 

refer to the contribution of transverse magnetic (TM) and transverse electric (Ti?) 
modes respectively. Here 



AT f (iO = V^-l + ei^O = p (5) 



where 



7i(*e) = y* 2 + ^(«) k 2 = ^(p 2 -l) (6) 

being the component of the field wave-vector parallel to the planar interface. The 
force per unit area is given by: 

The case of a free standing metal film is obtained with e 1 (c<j) = 62(00) = 1 and €3(0;) 
given by the Drude model. 

esH = 1 - ^. (8) 

where Q 3 is the plasma frequency of the film. One can easily verify that the use of 
the Drude model leads to a vanishing force in the two limits — > (infinitely diluite 
metal) and ^3 ^ 00 (ideal metal). Since the force cannot be identically zero for finite 
plasma frequency values, we expect that it has extrema in between these two limits. 

FigfS] displays the calculated force versus plasma frequency curves for a 50nm film. 
We have chosen uj t comparable with those obtained by fitting optical properties and 
from sample metal resistivity data [22]. For comparison we give in the same figure the 
curve obtained using the plasma model. It is seen that the general shape of the curve is 
the same: the force is attractive and shows a maximum of intensity (a minimum in the 
curve) and a long tail at high frequencies. However the inclusion of a finite relaxation 
frequency leads to significant modifications in the force intensity both at low and high 
frequency and causes a shift in the position of the maximum. 

In particular we notice that, on increasing the relaxation frequency, the force decrease at 
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Figure 2. Free standing film. Force on film boundaries as a function of film 
plasma frequency for different relaxation frequencies. lo t = continuous line, u T — 
lQ 1A rad/ sec long dashed line, uj t = 10 15 rad/ sec dot dashed line, uj t = 5 • 10 15 rad/ sec 
short dashed line. 




Figure 3. g factor of equation (fTT|) . gray line; / factor of equation (fT9|) . black line. 



low ^3 while it becomes significantly higher at high plasma frequency. The weakening of 
the force in the low plasma frequency range, can be understood using the small distance 
approximation to the Lifshitz formula, that can be applied when the film thickness d 
is much smaller than both the plasma wave length X p = 2txc/VL ? , and the relaxation 
wavelength A T = 2hc/oj t . In such conditions the force is due to TM modes only and is 
simply given by: 

h r (e 3 (z£) - l) 2 



F(d) 



(9) 



showing the typical d 3 behaviour of the van der Waals dispersion forces. Inserting the 
Drude dielectric function into this equation we get: 



where 



V2 



7T 
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ArcTan 
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x A 



(10) 
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and 

F « = -32^ < 12) 
is the force calculated for uj t = (plasma model) in the low thickness regime. 
Q s = Vtz/\/2 is the frequency of the free electron surface plasmon. Here we assumed 
uj t < \f2VL%, a condition that is commonly verified in real metals. The plot of the g 
function, presented in figj3] for a large interval of x values, shows that the correction 
to the force calculated in the ideal plasma model is always negative and gets more 
significant on increasing uj t . 

At high Q 3 values we do not have a simple formula like equation ffTUl) that allows 
to calculate the force analytically. The plasma wavelength is small and the condition 
d <C A p is not satisfied. However we can understand the reason of the increase in the force 
intensity compared with the plasma model using the following simple argument. The 
vanishing of the force as f2 3 goes to infinity as due to the exponential factor e - 2 £ K 3 d / c 
in Qtm and Qte (see equations (j3J) and (TjJ), since the fraction that multiplies the 
exponential goes to one both in the plasma and Drude model (this corresponds to the 
condition of reflectivity equal to one). We can write the exponential as e _27sd , for high 
f2 3 values we have 7! ~ £ 2 ^2+^ ? which shows that for the plasma model the exponential 
goes to zero faster that for finite u T values. This implies that both Qtm and Qte go to 
one more slowly as uj t increases, leading to a higher intensity of the force in the Drude 
model for high values. 

As previously observed [19] the occurrence of a maximum is a direct consequence of the 
retarded nature of the interaction. Indeed formula ([9]), which is obtained by neglecting 
retardation effects, does not give any maximum as a function of Q 3 . 
FigJH shows the calculated force per unit area in films of different thickness. The shape 
of the curves and the behaviour as a function of u T are basically the same of fig|2] and 
follow the same trend that has been reported previously for the plasma model [19]: the 
force decreases in intensity and the maximum shifts at lower frequencies on increasing 
d. Again we notice that the inclusion of the relaxation frequency can cause a significant 
enhancement of the force in the high f2 3 regime. Figj5] displays the plot of the force 
as a function of d for films of different plasma frequency. It is seen that the force, 
for relatively low plasma values (figj5] (a)) behaves as d~ 3 over a large thickness range. 
Significant deviations from this behaviour appear at very high relaxation frequencies 
only. At high plasma frequency (fig 15] (b)) different inverse power laws are present in 
the same interval of thicknesses, in agreement with previous findings [19] . 
The problem of the stability of a free standing film is commonly discussed in terms 
of competition between surface energy and surface stress and the volume strain energy 
[231 12^1 125] . Surface stresses can arise as a consequence of the different local environment 
of surface atoms with respect to those in the bulk: if they are tensile stresses they can 
originate a compressive stress within the film. If the compression is sufficiently large it 
can cause buckling of the surface with changes the lattice parameters of the outermost 
of atomic planes. Treating the film as a continuum one can perform the analysis of the 
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Figure 4. Free standing film. Force on film boundaries as a function of film plasma 
frequency for different relaxation frequencies. The notation and the oj t values are the 
same of fig'H] 



principal strains under the assumption that the film is in a state of biaxial stress [23J . 
For an isotropic system, taking the z-axis normal to the surface of the film, one has 
o~xx = o~ yy = a and a xy = a zx = a yz = o zz = and, neglecting quadratic terms, the 
equilibrium strains are given by: 

1 — 1/ v 

(-xx = £ yy = e || = a — e zz = e± = -2a — (13) 

where Y is the Young modulus and v the Poisson's ratio of the metal. In terms of the 
surface stress 

*W = 7 + ^ (14) 

the equilibrium condition requires 

= ~a (15) 

so that one can write the strains as 

£|l = -2.^ e± = ^ (16) 

here d is the film thickness and 7 is the surface energy per unit area. Typical calculated 
values of the surface stress are of the order of 1 -j- 5N/m [25J. Taking typical values for 
the elastic parameters one can estimate the strain values for d = lOnm to be of the order 
of 10~ 3 -T- 10~ 4 . To include the vacuum fluctuation forces, one has to carry out the same 
analysis adding to the surface and bulk elastic energy, the vacuum electromagnetic 
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Figure 5. Free standing film. Force as a function of the thickness for different 
plasma frequency and relaxation frequency values. uj t = 10 u rad/ sec diamonds, 
uj t = 10 15 rad/ sec crosses, ui T — 5 • 10 15 rad/ sec dots. Continuous line of figure (a) 
is ~ 1/d 3 " 6 , dashed line of figure (a) is ~ 1/d 3 , continuous line of figure (b) is ~ 1/d 4 
and dashed line of figure (b) is ~ 1/d 5 . 



energy of the film. This can be done along the same lines, leading to the following 
expressions of the strains 



u A-is F(d) /.x v F(d) , \ 

here F(d) is the vacuum electromagnetic force per unit area. As expected by a linear 
approximation, the strain induced by the surface stress and the vacuum electromagnetic 
force simply add to determine the strain condition of the system. Even if the 
contributions have opposite sign, the order of magnitude is different. For d = lOnm 
the maximum value of F(d), as shown in figJU is of the order of 10 A N/m 2 , leading to 
strains of the order of 10~ 5 -r- 10 -6 , which are too small to compete with the surface 
stress effects. Therefore, under normal conditions, the electromagnetic force does not 
contribute significantly to the free standing film morphology. For thicknesses of the 
order of few monolayers, the situation is less clear, since the electromagnetic force is 
expected to increase at least as <i~ 3 . However in such dielectric continuum theory 

is not adequate to represent the film optical properties, since size quantization effects 
become important. Their inclusion leads to a very different behaviour of the dielectric 
function and to an expression of the force where the derivative of the film dielectric 
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Figure 6. Film deposited on perfectly reflecting substrate. Force on film boundaries 
as a function of film plasma frequency for different relaxation frequencies. co T values are 
the same of ngf2] 



function with respect to d has to be considered. 

3. Deposition onto a perfectly reflecting substrate 

The situation can be significantly modified if one consider the case of deposited films. 
Previous work using the plasma model has shown that a metal deposited onto a perfectly 
reflecting substrate is subject to a positive pressure between the boundaries, that tends 
to increase the film thickness. This force turns out to be approximately double in 
intensity compared to the force on the free standing film [19]. 

Fig j6] reports the results of the calculation of the force acting on a 50nm film deposited 
on a perfectly reflecting substrate as a function of Q% for different uj t values. For 
comparison, the curve corresponding to the plasma model is also given. One can notice 
that the force is repulsive and its intensity is nearly doubled in comparison with the 
isolated film case. However the shape of the curves is basically the same: the force 
shows a maximum and a long tail at high frequencies. With respect to the plasma 
model, the inclusion of relaxation leads to an increase of the force in the high f2 3 region 
and a decrease at low frequencies. This behaviour, illustrated in figJEJ can be understood 
using the same arguments presented in the previous section. In particular, we notice 
that in the small d regime the force is simply given by: 





where 




(19) 



and 




-2F P1 



(20) 
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Figure 7. Vacuum energy second derivative plot for the Drude film deposited onto 
perfectly reflecting substrate. The notation and the uj t values are the same of figO 
The film thickness is d — 6nm. 

is the force calculated for the plasma model. The behaviour of the f(x) factor, 
reported in figJSl shows that the force is weakened compared to the plasma case. The 
strengthening of the force at high Q 3 is basically a consequence of the slower decay of 
the exponential in equations ([3]) and (j3J), as in the case of the free standing film. 
The d~ 3 dependence is characteristic of the small d regime only and it is appropriate for 
small plasma frequencies. As ^3 increase this approximation does not hold any more 
and, as in the case of the free standing film, different inverse power dependence are 
needed to represent the force behaviour. 

The presence of a positive pressure between the boundaries suggests that the force can 
be useful to stabilize the film. Indeed, under reasonable assumptions, one can show (see 
Appendix A) that a metal epitaxial film is stable against any small perturbation when: 



d 2 E(d) 



a 



1 



> TPT (21) 



dd 2 F 2 7 

where a is the mismatch lateral stress, Y is the Young modulus of the film, 7 is the 
surface energy and 9 is the second derivative of the vacuum fluctuation free energy 
with respect to the film thickness. 

The d value for which the two members of equation (1211) are equal gives the critical 
thickness d c , i.e. the maximum film thickness compatible with a flat surface morphology 
(see Appendix A). The critical thicknesses that have been reported in the study of 
metallic growth are of the order of few nanometers. In this condition and provided Q$ 
is not too large, one can use the small d approximation, which corresponds to the van 
der Waals force description. The energy second derivative turns out to be 

d 2 E 3V2hQ 3 



dd? 327Td 4 ' f \Q 3 ) ^ 
FigJ7]presents curves of d 2 E/dd 2 as a function of Q$ for a 6nm film at different relaxation 
frequencies in the small d limit. The plot allows to determine the range of f2 3 where a 
flat surface is stable, corresponding to the values of d 2 E/dd 2 above the horizontal line 
given by the second member of equation (l2"Tj) . which has been calculated taking typical 
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Figure 8. Critical thickness for a 6nm film deposited onto a perfectly reflecting 
substrate. The continuous and dot dashed lines represent respectively the plasma film 
and Drude films in the small d approximation, uj t = 5 • 10 15 rad/s. The long dashed 
and the dashed lines represent again the plasma and Drude films respectively, but they 
include also retardation effects. 



values of the parameters representing the elastic and surface properties of the film. We 
used the values o = 500MPa, 7 = U/m 2 , Y = 76GPa, reported in [10], changing 
these values simply shifts the horizontal line in the figure. The point where the second 
derivative curves cross the horizontal line is the minimum plasma frequency that allows 
for the flat configuration to be stable. It is seen that, given the film thickness, flat surface 
stability is possible above a threshold value. On passing from the plasma model, where 
d 2 E/dd 2 depends linearly on Q 3 , to the more realistic Drude description the stability 
threshold moves to higher Q% values. This indicates that, even adopting a simplified 
model of the substrate, the conditions of stability depend critically upon the description 
of the film optical properties. 

It is common practice to evaluate the vacuum energy contribution using the van der 
Waals energy. As we noticed before, this is a good approximation only when the film 
size d is much less that the plasma wavelength. If this condition is not verified, one 
obtains a rather inaccurate estimate of the critical thickness. To give an example we 
display in figJHl the d c values calculated with the full theory, and those obtained with 
the small d approximation. The comparison shows the inadequacy of the van der Waals 
description in the determination of film stability for a large interval of Q% values. Notice 
that both theories predict the existence of a finite d c at any Q% value. As shown in the 
next section, this result is peculiar and is a consequence of the perfect reflectivity of 
the substrate. Relaxation effects play a minor role and lead to differences in the critical 
values, more significant at large f^. 

4. Deposition onto a Drude metal substrate 

A more realistic description is achieved when the substrate dielectric function is 
described by a Drude model. Figj9] displays the calculated force per unit area for films 
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deposited onto a metallic substrate of fixed plasma frequency Q x as a function of fi 3 
for different combinations of relaxation frequency values. The thickness of the film is 
50nm. The curves are typical of the behaviour that can obtained from the calculations 
for this kind of systems. The behaviour is similar to the previously reported results 
on plasma model [19J. The force is repulsive in a limited range of frequencies and 
becomes negative when f2 x is smaller than the film plasma frequencies. This is strictly 
true when the relaxation frequencies of the substrate and the film are the same. If the 
substrate relaxation frequency is lower than the one of the film, the range of frequencies 
that corresponds to a repulsive force is larger and the force intensity increases, while 
in the opposite case the force curve crosses the horizontal axis at lower Q% values. To 
understand these findings we notice that the force calculated in the van der Waals, small 
d, regime, when the relaxation frequency is the same for the film and the substrate is 
given by: 

h n§(n*-nl)r 1 m . „ , u T 



— ArcTani 



cor 



1 (l-ArcTan{ , ^ 2 )) (23) 



10 



where Q s = Q 3 /\/2 and Q = a/ (fif + Q^)/2. Clearly this approximation predicts a 
repulsive force if ^3 < tt\ and an attractive one if ^3 > Qi in agreement with the result of 
the full calculation. This conclusion appears to be quite general and independent of the 
particular model of dielectric function adopted to describe substrate and film dielectric 
functions. It is consistent with the behaviour of the London dispersion force between 
dissimilar materials that has been reported in the literature, mainly with reference to 
ceramic or dielectric materials [27] . It has to be noticed that the force may change 
even by an order of magnitude when the relaxation frequency of the film is modified. 
This frequency is expected to be different in a film compared to the one measured in 
the bulk optical properties, due to the contribution of surface scattering and defects. 
This indicates that the use of the dielectric function of the bulk metal to represent the 
optical properties of the film may lead to serious errors in the determination of the 
vacuum force. Fig JTU] illustrates the behaviour of the force as a function of the substrate 
plasma frequency for a film with = 10 16 rad/ s, when the relaxation frequency is 
modified. The curves allow to understand how the force of the film is affected by 
gradually changing the properties of the substrate. These results like those of fig El are 
quite general since similar curves can be obtained by changing the film thickness and 
the other parameters. The examples reported in the figure show that the shape and 
intensity of the force can be significantly modified by varying the optical properties of 
film and/or substrate. It should be observed that in the real situation of film deposition, 
the optical properties of the film may be very sensitive to the experimental conditions: 
both and u T can vary significantly, well above the experimental errors, for films of 
the same metal prepared by different experimental techniques |20j . 
To discuss the stability problem we notice that, unlike the perfect metal substrate case 
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Figure 9. Film deposited onto a Drude substrate. Continuous line, w T 3 = ui T i = 
10 15 rad/ s; dashed line, uj T 3 = 10 15 rad/s and u T \ = 10 14 rad/s; dash dot lo t ^ = 
10 14 rad/s and ui t1 = 10 15 rad/s. 
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Figure 10. Film deposited onto a Drude substrate. The notation and the lo t values 
are the same of figE] 

illustrated in the previous section, deposition onto a real metal does not necessarily give 
a repulsive force. 

5. Stability in the van der Waals regime 

The issue of film stability is currently discussed using empirical forms of the interaction 
between the film boundaries based on the small d approximation [18J. To make reference 
to the previous work we first present the results obtained under this approximation, to 
be compared with those obtained including retardation effects that we give in the next 
section. To illustrate the importance of a realistic description of the optical properties 
we start discussing the outcomes of calculations using the plasma model i.e. neglecting 
relaxation effects both in the film and in the substrate. Fig{TT] illustrates the behaviour 
of d 2 E/dd 2 as a function of Q\ for films of different plasma frequency and 6nm thickness. 
For ^3 = 2 • 10 15 rad/s no stability region is achieved since the curve never crosses the 
horizontal line. This conclusion is obviously dependent on the parameters chosen to 
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Figure 11. Van der Waals energy second derivative plot. In the case of plasma model 
we used = 2 • 10 15 rad/s, continuous line and ^3 = 10 16 rad/ s, long dashed line. 
Short dashed line and dot dashed line have been obtained using the Drude model 
with f^3 = 10 le rad/s and iv T 3 — 10 15 rad/s, u) r \ — \Q 1A rad/s or uj t ^ — 10 14 rac?/s, 
uj T i = 10 15 rad/s respectively. The film thickness is d = 6nm. 



represent the elastic and surface properties. On the other hand for Q 3 = 10 16 rad/s we 
find a threshold value above which the flat film is stable. The difference is caused by the 
fact that for small Q3 values even a perfect metal does not provide a force repulsive 
enough to overcome the instability due to the elastic stress. Indeed for Q\ — > 00 
the second derivative of the energy for deposition on a perfect metal in the small d 
approximation is simply given by: 

d 2 E(d) = Wmh ( , 

dd 2 32vrd 4 1 ' 

if this quantity is less than one can never achieve a stability region. 
In the plasma model the only relevant parameters are the plasma frequencies of the film 
and the substrate. It is useful to represent the results in the form of a stability diagram 
in the Q3Q1 plane, like those given in figfl2l Each curve in the diagram is obtained 
at fixed film thickness as the locus of the couples Q 3 Qi which satisfy the equilibrium 
condition ( 12TI) . The plane is separated into two areas corresponding to flat surface and 
rough surface condition respectively. Notice that upon increasing the film thickness the 
domain of f^f^i values that allows for the existence of a flat surface gets narrower. In 
figjl3]we report the calculated critical thickness as a function of ^3 for a fixed substrate. 
Notice that, unlike the perfect metal case, there is a finite range of film plasma values for 
which the flat surface growth can occur. As expected this range reduces as fli decreases. 
The maximum value of the critical thickness is in correspondence with the maximum 
intensity of the repulsive force. 

Turning to the more realistic Drude model we notice in the first place that the values of 
the second derivative curve are systematically lower that those obtained in the plasma 
model, as clearly shown in figHU Upon increasing the relaxation frequency the crossing 
point shifts to higher Q 3 values. Since u T depends primarily on the quality of the 
film, this result suggests that the stability condition can be significantly modified by 
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Figure 12. Stability diagram for the plasma film deposited onto a plasma substrate, 
retardation effects are not included. Different curves for different thicknesses have been 
represented, above them the film is stable in its wavy morphology, below them the film 
is stable in its flat morphology. 
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Figure 13. Critical thickness. For the plasma model we used ill = I0 16 rad/ s, 
continuous line and fii = 5 • f 16 rad/s, long dashed line. Short dashed line represents 
the perfectly reflecting substrate f2i — > oo. For the Drude model we used a plasma 
frequency of Sli = 5 ■ 10 16 rad/s together with lu T 3 = w T i — 5 • 10 15 rad/s, dot dashed 
line. 
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Figure 14. Stability diagram for the Drudc film deposited onto a Drude substrate with 
the same relaxation frequency of the film, retardation effects are not included. Different 
curves for different relaxation frequencies have been represented, to T — 5 • 10 15 rad/s 
dashed line, u> T = fO 15 r ad /s dot dashed line and lu t = continuous line. 

the surface and defect scattering processes that take place in the film. The stability 
diagrams plotted in fig{TU for a 6nm films illustrate how the inclusion of relaxation 
processes can vary the range of plasma values allowing for a stable flat surface. We 
notice from fig {13] that the critical thickness values are reduced compared to the plasma 
model. 



6. Effects of the retardation in the evaluation of deposited film stability 

The plot of the second derivative curve calculated for the plasma model including 
retardation as a function of the substrate plasma frequency for a 6nm film, presented in 
fig US, shows that, unlike the van der Waals case, where stability is ensured for all the Qi 
values above the crossing with the horizontal line, the full theory predicts that the flat 
surface condition exists for a limited range of fl\ values only. This result puts limitations 
to stability arguments based on a van der Waals description of the vacuum force, since 
it shows that the small d approximation is not adequate even at nanometric thicknesses, 
if the substrate plasma frequency is high enough. The range is smaller when relaxation 
effects are included into the theory. The stability plots in the Q 3 Qi plane, given in figjl"6| 
show the modifications occurring in the range of film-substrate parameters compatible 
with the flat surface growth, when both retardation and Drude effects are included 
in the calculations. The reduction in the values of the critical thicknesses for a given 
Qi is illustrated in figJTTl It is clear from this analysis that the role of the vacuum 
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Figure 15. Vacuum energy second derivative plot including retardation effects. 
Continuous line represents the plasma model with = 10 16 rad/s. The Drude 
model, represented by the dashed line, is obtained with the same plasma frequency 
and ui T 3 = 10 15 rad/s and w r \ = 10 14 rad/s. The critical thickness is d = 6nm. 




Figure 16. Stability diagram with retardation effects. The dashed and gray 
continuous lines represent the plasma film deposited onto plasma substrate with and 
without retardation effects respectively. The dot dashed and black continuous lines 
represent the Drude film deposited onto Drude substrate with the same relaxation 
frequency lu t — 5 • 10 15 rad/s with and without retardation effects. 



forces in determining the stability if deposited films depends crucially upon the optical 
properties of both the film and the substrate. If a reasonably accurate description of 
the dielectric functions is not available, then it is hard to draw reliable conclusions on 
the expected surface morphology. Considerations based on simple empirical expressions 
of the force based on the van der Waals approximation may lead to overestimate the 
range of stability, if the conditions that make the small d approximation applicable are 
not met in the system under consideration. 
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Figure 17. Critical thickness including retardation effects. The gray lines represent 
the plasma and Drude film of fig ll3l in the small d limit. The continuous and 
dashed lines represent again the plasma and Drude models respectively, but including 
retardation effects. The substrate plasma frequency is fix = 5 • 10 16 rad/s. 

7. Conclusions 

The results we have presented illustrate the possibilities that metal films deposited on 
metal substrates be stabilized by vacuum fluctuation forces. We have shown that there 
is a variety of situations where these forces may play a role and we have expressed the 
constraints that have to be satisfied, in order to allow the maintenance of flat surface 
condition in the film growth, in terms of the parameters entering in the description of 
the optical properties of both the film and the substrate. 

The theory is based on a continuum description of the mismatch stress and a dielectric 
approach to the film and substrate optical properties. The application to a real system 
requires a detailed knowledge of the film optical properties. The requirement that the 
force be repulsive as well as its intensity may or may not be satisfied for a film of 
the same material depending upon the value of the relaxation frequency, the change in 
electron density that takes place in the pseudomorphic overlayer, the surface scattering, 
etc. A definite assessment on the role of the vacuum force on a specific system can be 
made only when these effects are properly accounted in the film dielectric function. 
The present theory can be improved along two main lines of development. One can 
improve the continuum approach (i) by adopting a more realistic description of the 
elastic energy in the study of deposited film stability with the inclusion of strain 
modifications in the substrate and surface stress effects [T7J [18], which have been 
neglected in the present paper; (ii) by including interband transition effects into the 
dielectric function, which are appropriate when transition metal or noble metal are 
considered [28]. One can also extend the theory to deal with the case of metal film onto 
semiconducting substrate. 

For very small thickness, i.e. ultra thin films with d < 1 4- 5nm, the theory has to 
be modified to include size quantization effects in the dielectric function [29]. This is 
expected to lead to significant modifications in the vacuum force, due to the change 
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in the film dielectric constant caused by the discretization of the electron energy 
levels. Moreover for ultra thin film the surface energy is expected to depend upon the 
thickness. Such modifications are consistent with the so called electronic growth model 
that has been frequently used to describe metal growth onto semiconductor substrate 



Appendix A. Macroscopic model for the film stability 

We consider the system sketched in fig. IA1I where a film of thickness d is deposited 
on an elastic substrate thick compared to the film. We assume that, under suitable 
conditions, the film can be grown as a single crystal on the substrate surface. If the film 
lattice constant is identical to the one of the substrate and the substrate-film interface 
is perfectly ordered, than the film will grow with zero macroscopic strain. On the other 
hand, in case of lattice mismatch, and still assuming a perfectly ordered interface, the 
film will grow continuining the substrate structure and will be strained in the plane 
of the interface, fig. IA1I (b). The film will tend to relax this elastic energy through 
dislocation motion or by changing the shape of the surface through mass transport, fig. 
IA1I (c). Here we report the analysis where also the second mechanism is present. It is 
generally accepted that at T = 0°K it is the competition between the elastic energy 
and the surface energy that determines the morphological stability of deposited film: 
indeed when a flat film surface is perturbed into a wave shape, as indicated in the fig. 
IA1I (a), elastic energy decreases, but surface energy increases. Surface roughening and 
island formation are expected to occur as a consequence of this competition, when the 
flat surface cannot be stable [121 II5[ 1121 1151 153] . 

However, when the size of the epitaxial film is nanometric, as in micro- nano-scale 
devices, other forces come into play. Long range dispersion forces, in particular the 
vacuum fluctuation force between the interface and the film free surface, have been 
considered as a possible source of stability in the system. The analysis is carried out for 
a sinusoidal surface considered as a perturbation of a reference perfectly planar surface, 
it is based on perturbation theory and on the assumption that the elastic constant of 
the film and the substrate are identical and in the absence of a stabilizing influence of 
the surface stress. For the two dimensional model sketched in fig JAll (a), where x is 
the surface position and z = is the reference plane and the sinusoidally wavy shape is 
described by a cosine function of amplitude q and wavelength A: 



one can show [34] that the stress concentration (stress for unit area and unit length 
normal to the surface) is related to the stress concentration a for the planar surface 
configuration by the equations 



[M EU [321 E3J. 




(A.l) 




(A.2) 
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Figure Al. Geometric parameters for the wavy film model and different possible 
stable configurations. 



tt^^fl-^).-^^) (A.4) 

which shows that at the surface the stress concentration is minimum in the peaks of 
the surface wave and maximum in the valleys and that in deeper locations away from 
the surface the stress is attenuated exponentially with a characteristic length A/27T. 
Assuming that d <C A one has, to the first order, that only the tangential stress 
contributes to the elastic energy 



a'(x) = a 



4vrg /2tt 
1 H — - — cos I — x 
X V A 



(A.5) 



and one can calculate the change in the elastic strain energy density on a wavelength as 

mi: 

AU el = -a 2 q\/Y (A.6) 

where Y is the Young modulus. Notice that AU e i does not depend upon the wavelength 
and, being negative, gives rise to film instability for any kind of harmonic perturbation 
of the surface. 

The change in the surface energy density can be determined by multiplying the surface 
tension 7 by the change in the surface length, which for small amplitude perturbations 
(q <C A) leads to the expression 

AU surf = 7 gV/A (A.7) 
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showing that the surface energy increases as the surface is perturbed harmonically. The 
comparison shows that, independently on the wave amplitude and the film thickness d, 
the flat film surface becomes unstable when: 

A > 7rr 7 /a 2 (A.8) 

the second member of this disequation giving a critical length below which the system 
is stable. Inclusion of the vacuum fluctuation force leads to an additional contribution 
in the energy density that can be calculated to second order in the surface modifications 
as: 



d 2 E(d) 



q 2 X 



d 4 



(A.9) 



a 4 hY 2 (A. 10) 



dd? 

Adding this term to the previous ones and imposing the equilibrium condition that the 
total energy change be zero, one arrives at the following conclusions: 

• the second derivative has to be positive, corresponding to a repulsive force on the 
film; 

• the film is stable provided its thickness d is lower than the critical thickness d c 
defined through the relation 

d 2 E(d) 
dd 2 

For film of very small size, one can assume a van der Waals expression of the vacuum 
energy 

AU vac = -Hq 2 \/8ird 4 (A. 11) 

where H is the Hamaker constant. One obtains the condition 

which shows that the film can be stable only for negative Hamaker constant (repulsive 
force) and it is inversely proportional to the stress concentration [10l [11] . 
Extension of the theory to three dimensions does not modify appreciably the stability 
condition is a single wave perturbation is considered. In particular equation flA.101) is 
still valid provided that one replaces Y with Y/(l — u 2 ), where v is the film Poisson's 
ratio [17] . 
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